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Abstract 

In models with universal extra dimensions (i.e. in which all Standard Model fields, 
including fermions, propagate into compact extra dimensions) momentum conserva- 
tion in the extra dimensions leads to the conservation of Kaluza-Klein (KK) number 
at each vertex. KK number is violated by loop effects because of the orbifold im- 
posed to reproduce the chiral Standard Model with zero modes, however, a KK 
parity remains at any order in perturbation theory which leads to the existence of a 
stable lightest KK particle (LKP). In addition, the degeneracy in the KK spectrum 
is lifted by radiative corrections so that all other KK particles eventually decay into 
the LKP. We investigate cases where the Standard Model lives in five or six dimen- 
sions with compactification radius of TeV -1 size and the LKP is the first massive 
state in the KK tower of either the photon or the neutrino. We derive the relic 
density of the LKP under a variety of assumptions about the spectrum of first tier 
KK modes. We find that both the KK photon and the KK neutrino, with masses 
at the TeV scale, may have appropriate annihilation cross sections to account for 
the dark matter, Qm ~ 0.3. 



1 Introduction 



One of the most exciting open questions on the interface between particle physics and 
cosmology is the nature of the dark matter. In fact, observations indicate that most of 
the matter in the universe is dark, and cosmological evidence has accumulated to provide 
independent confirmations that a large part of the Dark Matter (DM) is non-baryonic. 
Recent measurements of the Cosmic Microwave Background anisotropy combined with 
measurement of the Hubble parameter suggest a flat universe in which 30% of the energy 
density is due to non relativistic matter and only 4% is due to baryons, consistent with 
measurements from clusters and Big Bang Nucleosynthesis (see for recent reviews). 
In this paper we will use the value derived by Turner from combining all the current 



in which VLm is the matter density of the universe expressed as a fraction of the critical 
density for a flat universe, h is the normalized expansion rate (H = 100 h km s _1 Mpc _1 ). 

Observational evidence for DM has been building but we still have no solid clue as to 
its identity. Various candidates have been suggested and the theory of structure formation 
provides indirect evidence about some of its properties, strongly hinting that it is weakly 
interacting and non-relativistic at late times. In other words, it is cold dark matter 
(CDM). The standard model (SM) of particle interactions, while describing remarkably 
well the results of collider experiments, does not contain a suitable dark matter candidate, 
and thus it is necessary to consider extensions. There are essentially two well motivated 
DM candidates in this context: WIMPs and axions. WIMPs (Weakly interacting massive 
particles) were in thermal equilibrium with the Standard Model particles in the early 
universe. With masses in the 10-1000 GeV range and weak scale cross sections (cr ~ 10~ 9 
GeV -2 ) they would have fallen out of equilibrium such that their relic density today 
would correspond to n wiMP ~ Oil). Axions, originally postulated to address the strong 
CP problem, would not have been produced at thermal equilibrium (but through the 
decay of axionic strings or domain walls for instance). Their mass is constrained by 
astrophysical and cosmological arguments to lie in the range m ~ 10~ 5 -10 -2 eV. 

The most extensively studied DM candidate is the LSP (Lightest Supersymmetric 
Particle), a stable particle in supersymmetric (SUSY) models with conserved R-parity, 
which, in most SUSY scenarios, is the neutralino and is a typical WIMP. A broad range 
of experiments around the world are underway for detecting WIMPs, both through direct 
WIMP-nuclear scattering experiments and through indirect searches such as detection of 
cosmic flux from dark matter annihilation in the galactic center. Current searches are 
already exploring the parameter space of SUSY WIMPs. Unfortunately, SUSY models 
lack predictability. They contain a huge number of free parameters and one has to make 
several assumptions to reduce this number, for example by assuming a model to describe 
how supersymmetry is broken and how the effects of the supersymmetry breaking are 
communicated to the superpartners of the SM fields. For example, predictions for the 



data @: 




0.33 ±0.035, 
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cosmic flux from annihilation of the LSP in the center of the galaxy can vary over orders 
of magnitude when scanning SUSY parameter space. 

While the LSP is very well theoretically motivated, since the identity of the DM 
particles remains unconfirmed we should examine alternative possibilities. On the other 
hand, the DM issue sets important constraints on model building in particle theory. For 
any extension of the Standard model predicting the existence of a stable particle, one 
should compute its cosmological relic density to check whether it naturally accounts for 
DM or if it leads to overclosure of the universe in which case the model, or at least the 
cosmological picture associated with it, has to be revised. 

The issue when searching for a dark matter candidate is to find a stable particle. 
There are two options: 1) The particle essentially does not interact with the Standard 
Model particles, has a very small decay rate, and therefore is stable on cosmological 
scales. 2) The particle is coupled to the SM. In this case, there must be a symmetry 
to guarantee its stability. In the case of the LSP, there is R-parity to guarantee the 
stability. In this paper, we study a new DM candidate: the LKP (Lightest Kaiuza-Klein 
Particlef]) which interacts with SM particles and is stable because of a Kaluza-Klein 
parity. The LKP arises in a generic class of models in which all fields propagate in extra 
dimensions. The next section is devoted to explain these models, and the particle physics 
context. The model has many attractive features, including the fact that a relatively 
small number of parameters are sufficient to describe the LKP. Essentially one: its mass, 
which at tree level is the inverse of the compactification radius. In Section [5] we review the 
standard relic density computation. Our major work has been to calculate annihilation 
(and coannihilation) cross sections for the LKP in two cases. In the first case, the LKP 
is a Kaluza-Klein photon (Section |J), while in the second, it is a Kaluza-Klein neutrino 
(Section ||). We also study the effect of coannihilation in Section ||. Finally, Section [7] 
summarizes our results and discusses open questions which stimulate further work on the 
subject. Technical details are presented in the appendices. 



2 Universal Extra Dimensions 

Universal extra dimensions (UED) postulate that all of the SM fields may propagate 
in one or more compact extra dimensions |J. This is to be contrasted with both the 
brane world scenario || where the SM fields are constrained to live in three spatial 
dimensions while gravity can propagate in the bulk, and intermediate models |7j in which 
only gauge bosons and Higgs fields propagate in extra dimensions while fermions live 
at fixed points. However, there is significant phenomenological motivation to having 
fermions and gauge bosons living in the bulk, including motivation for three families from 
anomaly cancellation M, attractive dynamical electroweak symmetry breaking (EWSB) 
IPH p], ) (supersymmetric) models in which the Higgs mass is a calculable quantity 
||13|| , preventing rapid proton decay from non-renormalizable operators [|9,14], orbifold 



1 The first discussion on the relic density of stable Kaluza-Klein particles (referred to as "pyrgons") 
was made by Kolb and Slansky || . It was later alluded to in Q] . 
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breaking of the parity in left-right symmetric models \TE\ , and (through the mechanism of 



fermion localization) natural explanations for the observed fermion masses and mixings 
fLU , 1161 , 117, 18|,|19|. In this article we discover a new motivation for the UED scenario: to 
provide a viable dark matter candidate. 

The new feature of the UED scenario compared to the brane world is that since there is 
no brane to violate translation invariance along the extra dimensions, momentum is con- 
served at tree level leading to degenerate KK mode masses at each level and conservation 
of KK number in the interactions of the four dimensional effective theory. This statement 
is broken at the loop level, where the fact that the extra dimensions are compact leads to 
(calculable) violations of the full Lorentz symmetry [21], and as a result shifts the masses 
of the KK modes away from their tree level values. 

Further violations result by applying orbifold boundary conditions in order to remove 
unwanted fermionic degrees of freedom. These lead to loop contributions that are log 
divergent [22| in the effective theory, thus signalling that they cannot be computed but 
must instead be treated as inputs. They further correct the KK mode masses and break 
conservation of KK number to conservation of KK parity^, provided the terms induced at 
both of the orbifold fixed points are equal. Whether this will be true or not depends on 
the details of the compactification dynamics and the UV completion of the theory, but 
the assumption is self-consistent in the sense that if it is true at one scale, the cut-off scale 
for instance, it remains true at any scales since radiative corrections induce equal terms 
on both boundaries. The resulting theory has interactions only between even numbers of 
the odd-number KK modes. This conservation of KK parity implies that the lightest first 
level KK mode (LKP) cannot decay into SM zero modes and will be stable, in analogy 
with the lightest super-partner in a supersymmetric theory which conserves i?-parity. 
Thus, UED is the first extra dimensional scenario to predict a candidate particle for dark 
matter. A further consequence of KK parity, that KK modes must be pair-produced, 
leads to interesting collider phenomenology |||23 . 



At tree level the KK particles of a given level are predicted to be degenerate with 
masses n/R where R is the size of the compact dimension and n is the mode number. 
However, at loop level there are both calculable and incalculable corrections 121], GDI. We 



follow the perspective of Ref. [21] and treat the divergent corrections as perturbations 



on the 1/R masses of the KK modes. This assumption is self-consistent though not 
completely general, and could occur, for example, if for some reason the underlying theory 
causes them to vanish at the cut-off scale (A < 50i? _1 ). This prescription was employed 
in Ref. [|2T| and results in small (loop-suppressed) corrections to the KK mass spectrum 
induced by renormalization group evolution from A to 1/R. However, we do not strictly 
wed ourselves to the particular choice of the divergent corrections made in 21 ] , but instead 
allow ourselves the freedom to adjust these terms independently in the effective theory. 

For the LKP to be a well-motivated dark matter candidate, it should be electrically 
neutral and non-baryonic. Thus, the most promising candidates in the UED picture are 
first level KK modes of the neutral gauge bosons (analogues of the KK modes of the 

2 KK parity can be seen as the combination of a translation by irR with a flip of sign of all odd states 
in the KK Fourier decomposition of the bulk fields. 
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photon and Z), and the KK neutrino, z/ 1 '. One could also consider the first KK mode of 
the graviton, though this case seems less promising because its very weak gravitational 
interactions would imply that it will annihilate much less efficiently and could easily 
overdose the universe. Since similar incalculable loop corrections render the graviton mass 
a separate input of the theory, we may simply consider that the graviton is heavier than 
the LKP, such that at the time scales of interest to us all of the KK gravitons have already 
decayed into the LKP and zero modes. Alternately, one could consider a "deconstructed 
model" p3]] in which the extra dimension is represented by a chain of gauge groups and 



thus there need not be KK modes of the graviton. A simple two-site model can successfully 
reproduce the physics of the first level KK modes, and would be sufficient for our purposes. 
More convincingly, there is actually an argument for having the KK graviton heavier than 
the lightest KK neutral gauge boson. While the KK graviton receives negligible radiative 
corrections in contrast with the KK states of the SM particles and therefore has a mass 
equal to R~ l , one can easily see from Ref. j^TJ (see next paragraph) that the lightest 
state which diagonalizes the (B^\ W$ [ ) mass matrix is lighter than R~ x for R~ x > 800 
GeV. This is because receives negative radiative corrections. So even in the absence 
of boundary terms, the LKP, while being nearly degenerate with the KK graviton (the 
mass difference is less than 0.1 %) is indeed the KK photon for > 800 GeV. For 
smaller compactification scales, one can always argue that any tiny negative boundary 
term correction to the B^ mass will make the KK photon lighterfj. Our guideline is 
phenomenology. The only interesting and plausible cases would correspond to the LKP 
being a neutral weakly interacting particle. 

In the gauge boson sector, EWSB induces mixing between the gauge eigenstates, 
and in analogy with the familiar effect for the zero modes which produces the photon 
and Z boson. Including tree level contributions, EWSB effects, and radiative corrections 
to the masses, the mass matrix in the [B^ n \ W3 ) basis is [21] , 



- 1 2„,2 i 1„.„_„,2 \ 

(2) 



\gi92v 2 i + \gW + SM$ 



where g\ and g 2 are the U(l) and SU(2) gauge couplings, respectively, v ~ 174 GeV 
is Higgs vacuum expectation value (VEV), R is the radius of the extra dimension, and 
SMf and <5M| are the radiative corrections to the B^ and masses, including the 
boundary terms. In the absence of the radiative corrections, the mixing between the 
KK gauge bosons would be the same as that for the zero modes, and one would have KK 
modes of the photon and Z with the same Weinberg angle as the zero modes. The radiative 
corrections will generally disrupt this relationship, and each KK level will generally have 

3 One might be worried about the late decay of the KK graviton into the LKP. Indeed we know 
that unstable TeV relics may be dangerous if their lifetime exceeds 10 6 s because of their effects on the 
primordially synthesized abundances of light elements. However, this is true for large relic densities. In 
the case of (massless as well as KK mode) gravitons, we expect a very suppressed relic number density 
because of their very weak coupling. On the other hand, theoretical predictions of relic gravitons are 
necessarily subject to large uncertainties. A reliable estimate is difficult to obtain since it depends on the 
complicated dynamics of preheating and on the specific inflationary model considered. 
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two neutral bosons which are different mixtures of and As explained above, 

from an effective theory point of view 5M 2 and 5M| are separate inputs for the UED 
theory, but it is self-consistent to imagine that they are small and the resulting corrections 
to the tree-level n/R masses are modest. Within this framework one could also imagine 
that it is natural to expect 5M| > 5M 2 because g 2 > g\ and 5M 2 is further enhanced by 
larger group factors. For simplicity, we work in the limit 5M| — SM 2 3> gig2V 2 , so the 
mixing angle is effectively driven to zero by the large diagonal entries in the mass matrix^. 
Within this framework, one expects the LKP to be well-approximated as entirely 
It thus couples to all SM fermions (and the Higgs) proportionally to their hypercharges 
with coupling gi, and is approximately decoupled from the gauge bosons. 

Similar corrections apply to the Kaluza-Klein modes of the fermions, and generically 
their masses are also independent parameters of the theory. If one follows the prescription 
that the lightest particles are those which undergo only the U(l) hypercharge interaction, 
the lightest KK fermion would be the right-handed electron, . Note that the subscript 
U R" refers to the fact that it is a KK mode of the right-handed electron (and thus is 
an SU(2) singlet) as opposed to its chirality; it is a massive Dirac fermion with both 
right- and left-handed polarizations. In Section |], we consider the case in which is 
substantially heavier than B^ l \ and thus irrelevant in terms of its relic abundance while in 
Section |6.1| , we also consider the case in which is only slightly heavier than and 
thus coannihilation effects can be significant. If one relaxes the restriction that the fields 
which experience the SU(2) interaction are heavier than those which only experience the 
U(l), one could also consider the W^p or the as the LKP. We consider the case in 
which !/« is the LKP, including a variety of coannihilation channels, in Section y. 

We continue to consider the zero mode gauge bosons in terms of their well-known mass 
eigenstates, 7, Z and W ± , however we simplify our results by neglecting all EWSB effects, 
which correct our results at most by v 2 R 2 . Thus, we consider all of the SM fermions and 
gauge bosons as massless, and include the full content of the Higgs doublet (including the 
would-be Goldstone bosons) as massless physical degrees of freedom. This means that we 
neglect some processes, such as B^B^ — > zero mode gauge bosons all together, because 
they are v 2 R 2 suppressed compared to the dominant decay modes. It also means that we 
can choose to describe the neutral zero mode gauge bosons either in the Z 7 or in the 
W-P basis, as is convenient for the problem at hand. This approximation will be 
further motivated below, where we find that the favored regions of parameter space for 
dark matter have the mass of the LKP on the order of 1 TeV, much greater than v ~ 174 
GeV. 



3 Density of a Cold Relic Particle 

In this section, we review the standard calculation of the relic abundance of a particle 
species (see [[25|,^6| for more details) denoted Z which was at thermal equilibrium in the 

4 This is not so different from the situation in {H), for which sin 2 0$ < 0.01 for 1/R > 600 GeV. 
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early universe and decoupled when it was non relativistic. The evolution of its number 
density n in an expanding universe is governed by the Boltzmann equation: 

^ + 3Hn = -(av)(n 2 -n eq2 ), (3) 

where H = (Mp/3M Pl fl 2 is the expansion rate of the universe, n eq the number density at 
thermal equilibrium and (av) is the thermally averaged annihilation cross section times the 
relative velocity. We are eventually concerned by a massive cold dark matter candidate, 
for which the equilibrium density is given by the non relativistic limit: 



n 



^ g m 3, \-^, ( 4) 



where m is the mass of the particle species in question. The physics of equation (|]) is 
the following: At early times, when the temperature was higher than the mass of the 
particle, the number density was n eq oc T 3 , Z annihilated with its own anti-particle into 
lighter states and vice versa. As the temperature decreased below the mass, n dropped 
exponentially as indicated in @ and the annihilation rate T = n(av) dropped below H. 
The Z particles can no longer annihilate and their density per comoving volume remains 
fixed. The temperature at which the particle decouples from the thermal bath is denoted 
Tp {freeze- out temperature) and roughly corresponds to the time when V is of the same 
order as H. 

Equation ([3]) can be rewritten in terms of the variable Y = n/s, Y eq = n eq / s where s is 
the entropy s = 2ir 2 g*T 3 /45. g* counts the number of relativistic degrees of freedom. From 
the conservation of entropy per comoving volume (sa 3 =constant) we get n + 3Hn = sY 
so that 

sY = -(av}s 2 (Y 2 -Y eq2 ) . (5) 



We now introduce the variable: 



x = Y ■ ( 6 ) 



In a radiation dominated era, 



4tt 3 ^T 4 1 dx 



H = ^V- , t = — -»• — = Hx, (7) 

A5MI, 2H dt K J 



and (Ef) reads: 



dY _ (av) s 



Y 2 _ Y eq2 



dx Hx 

As is well known, (av) is well approximated by a non relativistic expansion (obtained by 
replacing the square of the energy in the center of mass frame by s = 4m 2 + m 2 v 2 ): 

(av) = a + b(v 2 ) + 0((v A )) w a + 6 b/x. (9) 

We finally rewrite our master equation @ in terms of the variable A = Y — Y eq : 

A' = -Y eq ' - f(x)A(2Y eq + A), (10) 
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where 

f(x) = ] j^-mM Pl (a + Q b/x) x~\ (11) 

A simple analytic solution can be obtained by studying this equation in two extreme 
regimes. At very early times when x « xf = m/Tp, A' << Y eql and A is given by: 

yeqi 

A = ~ f(x)(2Y eq + A) ' (12) 
At late times, A ~ Y » Y eq and A' >> Y eq ' leading to 

A- 2 A' = -/(*). (13) 
Integrating this equation between x F and oo and using the fact that A XF » A^: 

L F f{x)dx 

We arrive at: 



Y- 1 = \J^rM Pl m xp\a + 3b/x F ). (15) 
The contribution of the Z particle to the energy density of the universe is given by 

Qz = pz/pc, (16) 

where p c is the critical density corresponding to a flat universe, 

p c = ml M Pl /8n = 1.0539 x 10" 5 /i 2 GeV cm" 3 , (17) 

The value for h is given in equation (p. p z is simply given by pz = fn z n z = m^So^oo, 
So = 2889.2 cm -3 being the entropy today. Finally, the contribution to Q from a given 
non relativistic species of mass m z is: 

2 1.04 x 10 9 x F 1 

tozhr w — , 18 

m pi \fg* (a + ob/x F ) 

where g* is evaluated at the freeze-out temperature. For our cases of interest, we will 
have freeze-out temperatures in the region of 50 GeV, for which we use g* = 92. Note 
that the mass m z does not appear explicitly in this expression. Its effect is hidden in the 
coefficients a and b (of dimension GeV -2 ) as well as x F - Therefore, all we have to do is 
to compute the annihilation cross sections, expand them in the non relativistic limit and 
extract the coefficients a and b. We must also determine x F , the freeze-out temperature. 
The freeze-out temperature is defined by solving the equation 

A(x F ) = cY eq (x F ), (19) 
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using the expression for A(x) at early times, c is a constant of order one determined 
by matching the late-time and early-time solutions. It can be chosen empirically by 
comparing to a numerical integration of the Boltzmann equation. Equation ( |l9l) leads to 

g m Mpi(a + 6b/x F ) \ 

which is solved iteratively. The result does not depend dramatically on the precise value 
of c which we will take to have the usual value c = 1/2. 



%f = In c(c + 




3.1 Including Coannihilation 

As pointed out in |2lj||, the derivation presented above needs to be readdressed in the 



case of coannihilation. Here, we briefly summarize the approach presented in |26[ to be 
followed in this case. Such situation occurs when there are particles nearly degenerate 
with the relic Z but with masses slightly greater than rriz- These extra particles are nearly 
abundant as Z and if the mass difference is smaller or of the same order as the temperature 
when Z freezes out, they are thermally accessible and their annihilation will play a major 
role in determining the relic abundance of Z. Let us label Zi, i = 1, N, these particles 
nearly degenerate in mass. Z\ is the LKP, Z2 is the next LKP, etc. We also denote X, X' 
any zero mode (SM) particles. Reactions such as ZiZj «-> XX' change the Zi densities 
rij and determine their abundances. Since all -Zj>i which survive annihilation eventually 
decay into Z\, the relevant quantity is the total density of Zi particles n = 5Zi=i n ii an d 
the Boltzmann equation for n can be rewritten with accurate approximation as [BH]: 



where 



^ = _ 3iJn _( (Jcftl; )( n 2_ n e g2)) (21) 



N 

^« = £^^(! + A,) 3 / 2 (l + Aj ) 3/ V^ +A A (22) 

ij 9eS 



Oij is the cross section for the reaction ZiZj — ► XX', gi is the number of degrees of 
freedom of Zi, Aj = (m^ — m\)jm\, and, 

N 

g cS = J29i(l + \) 3/2 e- xA *. (23) 

i 

Equation is of the same form as (0) and can be solved using similar techniques. The 
formula for xp becomes, 

1 If , n , g^rn Mpi(a^ + Qb c «/x F )\ 
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where g has been replaced by g eB and a and b by a oB and b cB , the coefficients of the Taylor 
expansion for a aH . The relic abundance now reads: 

„ , o 1.04 x 10 9 GeV" 1 x F , . 

^ /2 M PI (/ + 3/6/2;/) 

with _ _ 

I a = xp / a off x~ 2 , and /& = 2x 2 7 / & off x~ 3 dx. (26) 

We now apply this formalism to our dark matter candidates. As explained in section II, 
we will consider both cases Z\ = flW and Z\ = We have computed the annihilation 
cross sections of the LKP into any zero mode (SM) particle. We begin by ignoring 
coannihilation and focus on the and candidates. Coannihilation effects will be 
considered in section |6|. 

4 as the LKP without Coannihilation 

We now analyze the case in which the LKP is B^ l \ and all other KK modes are consid- 



erably heavier (roughly 10% or more |[26||), so they do not play a significant role in the 
final relic density of the The relevant cross sections for pairs of to annihilate 

have final states into fermions or into Higgs bosons. In the limit in which EWSB effects 
are neglected, there are no channels into vector bosons. For simplicity, we neglect the 
mass splittings between the LKP and the higher states in the cross sections, expressions 
for which may be found in Appendix [A[ 

From these cross sections, we derive the coefficients in the thermal average discussed 
above, finding, 

_ Anal ( 2y F + ZYh) 

(X — 



9m 2 KK 
iral (2Y> + 3Y H ) 



iam KK 

where Yp = 95/18 and Yjj = 1/16. Numerically, for tukk — 1 TeV, the effective annihi- 
lation cross section is a ~ 0.6 pb, and the annihilation is 35% into quark pairs, 59% into 
charged lepton pairs, 4% into neutral leptons, and 2% into Higgs. For different masses, 
the cross section falls 8jS fnKK and the relative importance of the various final states stays 
approximately constant. As discussed above, these results allow us to determine xp, and 
we find that it is a very slowly varying function of tukk, decreasing from xp = 26 for 
ttikk = 200 GeV to xf = 24 for ttikk = 2 TeV. Therefore, it is essentially a and b which 
control the ttikk dependence of Q B (i) . 

In Figure [l] we present the prediction for flgco h 2 as a function of the KK mass for five 
dimensions. In five dimensions, an upper bound ttikk ^1-9 TeV is set from the universe 
overclosure condition and to account for the dark matter (Q = 0.33 ± 0.035), we find 
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Figure 1: Prediction for Q B (i)h 2 as a function of the KK mass (when neglecting coannihila- 
tion). The upper horizontal region delimits the values of Qh 2 above which the contribution 
from B^ to the energy density would overdose the universe. The lower horizontal band 
denotes the region Q = 0.33 ± 0.035 (using h = 0.69 ± 0.06) and defines the KK mass 
window if all the dark matter is to be accounted for by the B' 1 ' LKP. 

that the KK mass must lie in the range tukk ~ 900 — 1200 GeV, with a corresponding 
freeze-out temperature of order Tp ~ 36 — 48 GeV. These results are slightly above the 
experimental bounds on universal extra dimensions from precision electroweak data and 
collider searches (~ 350 GeV for one extra dimension ||), and imply that provided the 
fermion KK modes are not very much heavier than the future collider experiments 
will be able to study the region relevant for dark matter. 

5 z/ 1 ) without Coannihilation 

The situation is slightly more intricate in the case where the the neutrino is the LKP. To 
begin with, we now have a relic density composed of both the and its anti-particle, 
both of which annihilate among themselves as well as with each other. We assume that 
there is no cosmic asymmetry between particle and anti-particle in the analysis below. 
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If there were a large asymmetry generated before freeze-out, this effect could dominate 
the eventual relic abundance, and the computation below would have to be modified. 
We must also consider a larger number of annihilation processes, including final states 
of fermions, Higgs, and Z and W gauge bosons. The various channels are listed in 
Appendix || along with the necessary cross section formulae. 

We continue to consider the regime in which the other KK modes are considered light 
enough that we can neglect the mass splittings in the cross sections, but heavy enough that 
they do not result in a large modification of the final relic density. One would naturally 
expect the mass of the to be close to v^ , its weak partner. In fact any mass splitting 
between the two KK modes is an effect of EWSB, and could lead to dangerously large 



contributions to the T parameter p7[ . Such a contribution could be compensated by, i.e. 
a heavy Higgs boson f28| . As we will see below in Section |BT2] , including a degenerate 



will not substantially alter our results. 



5.1 One Flavor 

Our relic density is both v^ and (n z = n v (i) + with g eS 
annihilation cross section appearing in the Boltzmann equation is 



4) so that the 



B 



with 



vv) + a(v^v {1) 



W + 2 a(v {L, v 



X) 



(28) 



aiv^v^ 



+ a(v^v^ 



qq) + 

-> ZZ) + cr(u 



W) + aiv^v^ 



(i)p(i) 



w + w- 



i+r) 



(29) 



where the cross section into quarks contains a sum over all quark flavors, the cross section 
into neutrinos contains a sum into both the neutrino zero mode of v^ and the other 
flavors, and the cross section into charged leptons includes both the zero mode charged 
partner of v^ l \ and also the other flavors. Note that the matrix elements for annihilation 
into other flavors are different from those into zero modes of the same flavor. 

Proceeding as before, we expand the effective cross section in powers ofl/xp, obtain- 



ing, 



a 2 n (2724, - 281s^ + 154) 



a 2 7t (2744, - 1274, + 4 ) 



15364 



W C W m KK 



(30) 



where the vv portion of the result sums over all allowed final states, including 3 up- and 
down-type quarks, 3 charged and neutral leptons, ZZ and W + W~ weak bosons, and the 
Higgs doublet. For ttikk = 1 TeV, we have the effective cross section a cB = 1.3 pb, slightly 
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0.5 1 1.5 2 2.5 3 



Figure 2: Prediction for Vt u (\)h 2 as a function of the KK mass. The solid lines are for 
u^ 1 ' alone (in the one and three family cases) and the dotted ones correspond to the cases 
where coannihilation with degenerate is included. 

higher than the B^ case. a cH is composed 18% of the process u^u^ — > vv^ with the 
remaining 82% coming from — > X. This second contribution is roughly 41% into 
quarks, 7% / 9% into neutral/charged gauge bosons, 2% into Higgs, and 33% / 8% into 
charged/neutral leptons. 

Deriving the freeze-out temperature, we find that xf varies from 27 for itlkk = 0.2 
TeV to 25 for vtlkk = 2 TeV. z/ 1 ) therefore freezes out somewhat later than and 
thus has a smaller relic density. The higher effective annihilation cross section translates 
into a different prediction for the KK mass to account for the dark matter energy density: 
n^KK ~ 1-3 — 1.8 TeV and the overclosure limit is pushed up to 2.7 TeV (see figure 
0). Again, these values are within the reach of planned experiments such as the LHC, 
provided the colored KK mode masses are not significantly different from the mass of the 
KK neutrino. The value for xf is not much different from the B^ case, however, given 
the different m KK window, the freeze-out temperature is higher, Tp ~ 50 — 70 GeV. 
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5.2 Three Flavors 



If we consider three degenerate flavors of KK neutrino, the relic density is computed 
as the sum of the densities of all three species plus the sum of the corresponding anti- 
particles. Thus we have g cB = 12. In this case, the effective cross section contains three 
separate contributions identical to that considered above, for each species to annihilate 
among itself, and also additional cross-flavor channels such as v^up — > V1V2, Vx^V^} — > 
v{T>2-, etc through t-channel Z^ 1 ' exchange. There are also cross-flavor transitions into 
charged leptons through t-channel exchange. The relevant formulae may be found 
in Appendix [FJ. We continue to assume no cosmological asymmetries between particles 
and anti-particles, and further consider the case where there are no asymmetries between 
different flavors. 

The effective cross section becomes, 



a cB - — 



+ 2a(u[ 1] ui 1} -> v x v 2 ) + 2a(F5 1) F^ 1) -> V X V 2 ) + Aa^V^ -> X)\ (31) 

where we have assumed that cross sections for all flavors (and combinations of flavors) 
are equal. The cross-flavor annihilation channels are not as efficient as the same-flavor 
channels (about twice in size), and Xp is about the same as in the single flavor case. 
Thus, the net result is a larger predicted relic abundance for the same mass, as shown in 
Figure §. Thus, the region relevant to explain measurements is lower, ttlkk ~ 950 — 1250 
GeV, and the overclosure condition requires uikk ^,1-9 TeV. The freeze-out temperature 
in the relevant region ranges from 36 — 47 GeV. 



6 Coannihilation Results 

Coannihilation is expected to play a significant role when there are extra degrees of 
freedom with masses nearly degenerate with the relic particle. Experience with the su- 
persymmetric standard model indicates that large effects are to be expected when the 
heavier particles have masses within about 5% of the LSP. The radiative corrections to 



the KK spectrum under the prescription of Ref. [51] indicate that quark and gluon KK 
masses can be shifted by twenty percents. Weak gauge bosons also receive corrections (at 
tree level) larger than five percents so that the only particles which will be considered 
as nearly degenerate with the LKP are the leptons. We will simply our analysis by con- 
sidering all higher Kaluza-Klein modes relevant for coannihilation to be degenerate, and 
leave the splitting between the LKP and next lightest Kaluza-Klein particle (NLKP) as 
an adjustable parameter. 

As motivated in Section ^, we will compute coannihilation channels in the two following 
situations: 

• is the LKP and is the NLKP with all other KK modes heavy enough that 
they do not contribute to coannihilation, 
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• is the LKP and is the NLKP (almost degenerate). 

The relative mass difference between the LKP and the second LKP is denoted by 
A = (m NLKP - m LKP )/m LKP . 

6.1 B^> Coannihilation with e$ 

In the first case, we consider as the LKP and as the NLKP, assuming no net 
asymmetry between the number of and e^. We consider both the case with one 
family of , and also the case of three degenerate families. With one family, the formula 
( p3|) for the effective number of degrees of freedom becomes, 

g eS = 3 + 4(l + A) 3/2 exp[-xA], (32) 

where we have used gj = 2,g B i = 3 and summed over B^\ v^ 1 ' , and Ta 1 \ The effective 
annihilation cross section is, 

fa* = gla(BWBW)+Ag B ,g f [l + Af 2 eM-xAHBWe^) 

+2g}[l + A] 3 exp[-2Ax] (*(e$e$) , (33) 

where we have assumed that the cross sections for annihilation of (B^ejJ ,5^e^) and 

( e R e R j £r £r) m t° zero modes are equal. 

The cross section a[B^ v 'B^ 1 >) is as derived before in Section [|. The cross section 

a{B^ep) proceeds into final states with zero modes of and eZ, and in the limit in 
which the Z mass is neglected can be equivalently described as a single process 
eB(°\ No vW~ final state occurs because e^, being a weak singlet, does not couple to 
the SU(2) bosons, and we neglect the tiny electron Yukawa coupling which would result 
in a e$° final state. There are also channels which convert e^e^ into fermions and Higgs 
through an s-channel B^ and into two B^ J s (or equivalently, into ZZ 77 and Z7 final 
states); and channels in which e^e^ exchanges a B« to become lepton zero modes. 
If more than one flavor of has a mass close to B^\ one also has channels in which 
different flavors of exchange a t-channel B^ and thus scatter into their corresponding 
zero modes. All of the needed cross sections are given in Appendix |Q. 

Our result when including almost degenerate with (A = 1%) is a higher 
LKP relic density than in the case without . Indeed, the self annihilation cross section 
of is not much higher than the one for B^ and the coannihilation cross section is 
significantly smaller (there are only two coannihilation channels while B^ and can self 
annihilate into all zero mode fermions). This situation is to be contrasted with the SUSY 
case where coannihilation between the neutralino and sfermions can be very efficient and 
significantly reduce the relic density. Here, we have more relics (both B^ and ) which 
essentially decoupled at the same time (and at roughly the same freeze-out temperature 
as was the case for B^ alone) and eventually the left over decay into B^\ This 
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Figure 3: Prediction for Vt B (i)h 2 as in Figure |l]. The solid line is the case for B^ alone, 
and the dashed and dotted lines correspond to the case in which there are one (three) 
flavors of nearly degenerate e^. For each case, the black curves (upper of each pair) 
denote the case A = 0.01 and the red curves (lower of each pair) A = 0.05. 

translates into a KK mass window slightly below the window obtained for alone. In 
Figure |3] we present the resulting relic abundance of B^ including both the one flavor 

— ( 1) 

and three flavors of , for two choices of A corresponding to 1% and 5% mass splittings. 
The curves become approximately degenerate with the B' 1 ' without coannihilation case 
when A > 0.1. In each case, the resulting tukk window shifts slightly downward because 
of the increase in the predicted relic density, favoring values between 600 — 1050 GeV, 
depending on the number of light flavors and the mass splitting. 

6.2 Coannihilation with 

As mentioned in the introduction of section 5, one should include in the calculation 
of the LKP relic density when assuming that the LKP is . Indeed, and are 
expected to be nearly degenerate, with tree level mass splittings on the order of the mass 
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of the charged lepton. For one family of leptons we have, 

g cB = 4 + 4(1 + A) 3/2 exp[-xA], 
and the effective annihilation cross section is: 



(34) 



g) ^(z/V 1 ) -> vv) + 2a(u^V^ 



+ 
+ 



2a{e 
4a(v 



e e 



X) 

->X)] (1 + A)V 2Ax 
+ X)] (l + Af/V lA } 



(35) 



where, 



tr^pM - X) 



+ cue 



aU 1 ) 



ZZ,Z 7 , 77 ) + cr(e^ ij e^ j -> W+W") 



(1)7,(1) 



i + r) 



+ cr(e 



(i)^(i) 



L e L 



qq') 



+ cr eVV 



(1)77(1) 



+ <r[e\'v 



(1)77(1) 



(1)77(1) 



W+7) 



+ cr (e}/u 



(1)77(1) 



For the three family case, we also include the cross flavor annihilation channels, 

g cS = 12 + i2(l + A) 3/2 exp[-xA], 



(36) 

(37) 
(38) 



and, 



„ 2 ~ 



3gj {2cr(z/( 1 )z/(D - w) + 2<t(iW 1 > - X) 

Z/ 1 Z/ 2 )+4cr(z/f ) z7( 1) ^X) 
r r j jlrr) r, r , — > X) 

,w..(i) > , A-fjn.s 1 ) > 



+ Aa{^^ 



+ 



+ 



+ M#VE ; - e-Ai") + ^(^V^ - X) 

Aaiefv^ -> X) 

^/i-I7)+4 ( T( Ai i 1) Z/( 1 )^X)l(l 



;i + a)V 



-2Aa; 



4<x(z/ 



+ 4<r(/x 



(i).(i) 



A) 3/2 e -,A| ; (39) 



with, 



,(1)^(1) 



e /i 



,(1)^(1) . 



v e v. 



(40) 
(41) 



Cr(/X?Z/(l) -> X) = (T^U/i-l/) + CX^V 

All relevant cross sections are given in Appendix [Cj From these results, we derive the 
freeze-out temperature and final relic density as a function of the mass splitting, A = 
(m (i) — m v {i) ) / m v (i) . We find that the result when A = is a small modification of that 

e L 

with 2/(1) alone. The modification of the freeze-out temperature is less than 1% over the 
relevant mass range, and the final relic density, shown in Figure |2|, is slightly decreased 
in the one family case and increased in the three family case. As A increases, the results 
rapidly return to the corresponding z/ 1 ) results. 
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7 Summary and open questions 



In this work, we have computed the contribution to the energy density of the Universe 
coming from the relic density of the Lightest Kaluza-Klein Particle in two cases: 1) The 
LKP is a Kaluza-Klein photon, 2) the LKP is a Kaluza-Klein neutrino. In models with 
universal extra dimensions (UED), such particle is stable and provides an interesting Dark 
Matter candidate whose mass is the inverse of the compactification radiusf] R. One could 
question under which general conditions is the KK parity preserved. There can indeed 
be violations of KK parity because of the states and interactions which are localized 
at the two orbifold fixed points. However, if the boundary Lagrangians are symmetric 
under interchange of the two boundaries, KK parity will be preserved and unaltered by 
radiative corrections. The boundary Lagrangian describing any coupling between states 
localized at the boundaries and bulk fields is invariant under KK parity when fields and 
couplings are the same at the two fixed points. In a string theory context, couplings 
between bulk modes and twisted sectors can be predicted. While the symmetry of the 
boundary Lagrangians is not generic, it is not either disfavoured and it is not difficult to 
imagine string theories which could preserve KK parity and still have all of the features 
necessary. The assumption that KK parity is preserved does not conflict decisively with 
any feature that a low energy theory derived from string theory must possess. 

We have assumed for simplicity that the Standard Model lives 4+1 dimensions. One 
could generalize our results to a higher number of extra dimensions^. The generalization 
to the six- dimensional case is reasonably straightforward. Let us consider for simplicity 
two extra dimensions with the topology of a torus, and equal radii. We impose either 
a Z 2 or Z^ orbifold symmetry, so our spaces are T 2 / Z 2 and T 2 / 'Z 4 . For T 2 /Z 2 our 
particles have now two KK numbers and there are two LKP degenerate in mass: B^ 1 ' ^ 
and B^ ' 1 ^ and two heavier stable states, B^ 1 ' 1 ^ and B^ 1 '^ 1 ^ whose tree-level masses are 



about 40% larger. Following the formalism developed in Section J7L, we see that for such 
large mass splittings, we expect the contribution of the heavier stable states to Qm to be 
exponentially suppressed compared to the light states, and thus we neglect them in the 
analysis. Because of KK number conservation (mod 2), at tree level the two LKP do not 
talk to each other and annihilate independently. The number of zero modes remains the 
same, therefore the relic density is just twice the one computed in the 5D case. The KK 
mass window to account for the DM is shifted to 650-850 GeV and the limit for overclosure 
is 1.3 TeV. For T 2 / Z± there are two stable particles with KK numbers (0, 1) and (1, 1). 
Again, the heavier state has a tree level mass about 40% heavier than the LKP, and 
we thus conclude that the prediction for Q M is approximately the same as the 5d case. 
The generalization to seven or more dimensions is more subtle because the number of 



5 Note that in models with TeV compactification scale, the higher dimensional Planck scale is ~ 10 14 
GeV for one extra dimension and ~ 10 11 GeV for two extra dimensions. 

6 One should also keep in mind that the computation of the relic density is a sensitive function of the 
expansion rate of the universe. The standard calculation assumes that the expansion rate is given by the 
conventional Friedmann equation evaluated in a radiation-dominated era corresponding to universe made 
of a gas of relativistic SM particles. Any deviation from this assumption will affect the relic density. 
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fermionic degrees of freedom is modified. One advantage of this model is that the physics 
is dominated by a single parameter: the size of the extra dimension, R. Interestingly, 
for the UED model to explain the Dark Matter with the LKP, we find that R typically 
has to be of the TeV scale, which is phenomenologically interesting and relevant at future 
colliders. 

Having checked that the prediction for Q, M is of the right order, the next step is to 
consider detection. Similarly to other WIMPs, the direct search for the LKP relies on the 
deposition of ~ keV recoil energy when the WIMP scatters from a nucleus in a detector. 
To study in more details the constraints on the LKP as the dark matter, the computation 
of the corresponding elastic scattering cross section between (or v^) and a nucleus is 
needed. This task is beyond the aim of this paper. In addition, indirect WIMP searches 
rely on the detection of 7 rays, charged particles or neutrinos from WIMP annihilation. 
There are two places where annihilation can take place: 

• In the Sun where the LKP may be captured and annihilation greatly enhanced. 
This will generate a neutrino spectrum. The prediction essentially depends on the 
competition between the gravitational capture of the LKP by the Sun and the LKP 
annihilation so we would need to know the details of the capture rates and of the 
propagation of the neutrinos from the core to the surface of the Sun to make any 
statement. 

• In the core of the Milky Way. LKP annihilation is important in the galactic center 
where the matter density is higher. To compute the resulting spectrum, one needs 
to know the reprocessing of the direct products of LKP annihilation. 

Among secondary products of annihilation in the galactic center are high energy 7 orig- 
inating via neutral pion decays (pions result from the hadronization of the directly pro- 
duced quarks) and the synchrotron radiation of e + e"~ pairs originating from the decays of 
charged pions in the galactic magnetic field. This requires the implementation of fragmen- 
tation functions. However, the flux of neutrinos coming from the direct LKP annihilation 
in the galactic center can be determined reasonably model-independently. These neutri- 
nos will not be reprocessed during their journey between the galactic center and us. To 
do that, we use a Navarro-Frenk- White profile for the Milky Way of the form: 



Po = 0.3 GeV cm" 3 is the local halo density, R$ = 8 kpc the distance between the Sun 
and the galactic center, a = 20 kpc some length scale. The LKP annihilation rate is then 
given by 



which leads to a flux of neutrinos at one TeV (assuming a mass of 1 TeV for the LKP) of 
4.4 x 10~ 12 cm" 2 s" 1 if the LKP is a KK photon, 2.3 x 10~ 10 cm" 2 s" 1 if the LKP is one 
flavor of KK neutrino and 2.8 x 10~ 9 cm -2 s _1 if the LKP is three flavors of KK neutrino. 




(42) 




(43) 
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Initial state 



Final state 
(zero modes of SM fields) 



Feynman diagrams 



fid) 



0* 




Table 1: Feynman diagrams for which we calculate the annihilation cross section of a KK 
photon into SM particles. s(x), t(x) and u(x) denote a tree-level Feynman diagram in 
which particle x is exchanged in the s-, t- and w-channel respectively. "Contact term" 
represents the scalar-gauge boson four point interaction. / denotes any zero mode fermion 
and is the scalar Higgs doublet. 

This is within the sensitivity of some future km 3 neutrino telescopes (like IceCube which 
will reach a sensitivity of about 10~ 12 in these units at a TeV). If the photon fluxes are 
of the same order in first approximation, experiments such as MAGIC which will reach 
about 10~ 12 at a TeV could set constraints on the mass of the LKP. While a large number 
of experiments are underway, this issue of Kaluza-Klein dark matter detection will be 
worth pursuing. 
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In this appendix, we summarize the annihilation cross sections for B^> into SM fields. 
In discussing the various processes it is worthwhile to remember some general features of 
the couplings of the first level KK modes we are considering. First, KK parity insures 
that only vertices with even numbers of first level KK modes exist. The two important 
types are couplings of two KK matter fields (fermions or Higgs boson) to a single zero 
mode gauge boson and coupling of a first level gauge boson to a KK matter field and a 
zero mode matter field. Recall that the left- and right-handed components of the zero 
mode fermions each have separate massive KK modes, whose couplings to the zero mode 
gauge bosons are vector-like. When a KK mode gauge boson couples to a KK fermion 
and a zero- mode fermion, there are generally projectors which insure that the zero mode 
fermion has the same chirality as the KK fermion. 



A Annihilation Cross Sections 
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In the limit in which electroweak symmetry breaking effects are neglected, pairs of 
can annihilate into either zero mode fermions, //, or pairs of Higgs bosons, 00*. 
Figures |] and | show the relevant Feynman diagrams. We approximate all first level KK 
masses (m) as equal, and ignore all "zero mode" fermion, scalar and gauge boson masses. 
The cross section for B^B^ — ► // receives contributions in which both the fj/ and 
fiP (in the case of neutrinos only u^') are exchanged in both the t- and u- channels. 
Summing/averaging over final/initial spins and integrating over the phase space of the 
//, the result may be written, 

a{B B -//) = _ — (44) 

where /3 is defined as, 




/3 = f-— (45) 

and for each fermion, ql = Q\Yl and Qr = giYjt. Note that the massless fermions in the 
final state have prevented interference between the graphs in which is exchanged and 
fx is exchanged. The factor N c sums over the different color combinations allowed in 
the final state, N c = 3 for quarks and N c = 1 for leptons. Thus, the sum over all three 
families of SM fermions results in, 

Q5 

NM + 9r) - 3gKYe\+Yl + Y* L +3(Y u \+Y: R + Yl+Y* R )) = -gl (46) 

Annihilation into Higgs, with Feynman diagrams shown in Fig. [5], proceeds through t- 
and it- channel exchange of <f>^ as well as the four-point interaction of B^B^<f><f>*. The 
cross section is, 

a(Z^5« ^ 00*) = 9 ^±- (47) 

where Y^ = 1/2 and a factor of two is included from the sum over the two complex 
fields in the scalar doublet. Note that this result implicitly includes the decay (after 
EWSB effects are properly taken into account) into longitudinal W and Z zero modes as 
well as the Higgs particle h, as all of these degrees of freedom are included in the scalar 
doublet. Neglecting EWSB there are no tree-level decays into gauge bosons; such decays 
are induced by EWSB, but are suppressed by powers of v 2 /m 2 and thus we neglect them 
here. 



B z/ 1 ) Annihilation Cross Sections 

The properties of the KK mode of the (left-handed) neutrino are assumed to be ap- 
proximately independent of the neutrino species. For simplicity, we do not consider the 
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Figure 4: Feynman diagrams for B^B^ annihilation into fermions. 

B 1 vWM/Vf|- - - § B 1 ^ ^p-> 4> B 1 ^ x (|) 

II 

II 

B'ww#-<--t* b 1 r ^--<-- f B 1 ^ "n ()>■ 




\ . * 



Figure 5: Feynman diagrams for B^B^ annihilation into Higgs scalar bosons. 



possibility of a sterile neutrino or its KK modes. For the purposes of this discussion, 
we assume a neutrino which is the weak partner of the left-handed electron; the results 
for the weak partners of the muon or tau are simply obtained by appropriately replacing 
the exchanged particles in specific processes. We continue to neglect fermion and boson 
masses, and ignore fermion mixing. 

The can annihilate with z/W into quark (and other family lepton) zero modes 
through an s-channel Z zero mode (Figure |6]). The cross section is given by, 



ff) 



N c g 2 z (9 2 l + 9r) (s + 2m 2 



where, 



9z 



(46 



(49) 



are the couplings of the Z° to u^'u^ 1 ' and 9l{r) are the standard zero mode couplings 
between the Z° and //, 



9l(r) 



SyyCw 



T 3 -Q f s 



(50) 



where T 3 is the third component of weak iso-spin of / and Q f its charge. N c accounts for 
the sum over final state color configurations, as before. 

Annihilation into zero modes of the charged lepton partner e + e~ proceeds either 
through an s-channel Z or a t-channel (Figure |7|), or into its own zero modes (uV) 
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Initial state Final state Feynman diagrams 
(zero modes of SM fields) 

i/WpM qq s(Z) 

0* s(Z) 

vv s(Z),t{B^,wP) 

e~ e+ s(Z), t(W + {1) ) 

Z° Z° t(vW), u(uM) 

W~W + t(e ( p), s(Z) 



V 



v x v 2 t(ZW) 



Table 2: Same as Table [I] but for annihilation of the KK neutrino. 




Figure 6: Feynman diagrams for annihilation into quarks or leptons of other 

families. 



through an s-channel Z zero mode or by exchanging a t-channel W 3 or (Figure |?p. 
In the limit in which we ignore the mass splitting between and B^\ the exchange of 
both can be summed into a single Z^ exchange. The cross section into either zero mode 
neutrinos of the same flavor or their charged lepton weak partners is, 

,2 _ L -^ f„ + 2 m?) 



11) 



9 Z glg L Ws + 2{2s + 3 m 2 )L] + g\{g\ 



t R ) 



2 „2 



32 71 P' 

gi [P (4s + 9m 2 ) + 8m 2 L] 
64 7r m 2 (3 2 s 



where g L = g z , 



log 



1-/9 ' 
l+P 



(51) 



(52) 



and g L ( R ^ are given as before. For annihilation into charged leptons, the g L m\ should be 

replaced by the charged lepton values from Eq. |50] and cjl now corresponds to the W± 
coupling to and e°, 

e 



9l 



V2s 



(53) 
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I 1 




I 1 



V 



Figure 7: Feynman diagrams for „(%(!) annihilation into zero mode leptons of the same 
family. 



Figure 8: Feynman diagrams for annihilation into scalar Higgs bosons. 

First modes of v^V^ can annihilate into zero modes of Higgs bosons through an 
s-channel Z zero mode (Figure ||) with cross section, 

48 Ti p s 2 

where the gz coupling of to a Z zero mode is as before, and there are two couplings 
of 00* to Z° from the charged and neutral entries of the doublet, respectively, 



9<f> 



SwCw 



(T 3 - Q, 



s w ) ■ 



(55) 



where the first (upper entry in the doublet) Higgs has charge = 1 and the second has 
Q - 0. 

Annihilation into ZZ are mediated by t- and w-channel , (Figure |9]) and has cross 
section, 

g| (2 [s 2 + 4m 2 s - 8m 4 ] ArcTanh[/3] - + 4m 2 ]) 
^ V " 87^ [bb) 

where gz is the coupling to the zero mode Z boson defined above. 

Annihilation into W + W~ zero modes includes t-channel exchange and s-channel 
annihilation through a virtual Z zero mode (Figure |T0|) . The cross section is, 



^(2/1)77(1) ^ W + W~) 



-5g 2 zww9z[ s + 2rr f\_ + 9w9z ww9z[l3s - 2m 2 L] 



2Att(3s 2 8ir(3 2 s 2 
g^[f3(s + 4m 2 ) + (s + 2m 2 )L} 
8tiP 2 s 2 



(57) 
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Figure 9: Feynman diagrams for 

^(i) F (i) 

annihilation into two neutral vector bosons VV . 
I 1 \\ W + I 1 > i pAAAAA^ W + 



W" / < W wwv W 



Figure 10: Feynman diagrams for annihilation into two charged vector bosons 

W+W~. 



where gz are the couplings to the Z zero mode as before, gw is the (vector-like) coupling 



between the W zero mode to and v^ x \ 



9w 



y/2sw 



and gzww is the Z-W + -W coupling between zero modes, 

c w 



9zww 



s w 



(5f 



(59) 



Furthermore, 

„(!)„(!) (77(1)77(1)) can annihilate into vv [yv] through t- and u-channel 
exchange of or (Figure |Tl|). The cross section for annihilation of two neutrinos 
pWi/i 1 ) — » W i s given by, 

gt(f3s(2s - m 2 ) + 2m 2 (4s - 5m 2 )ArcTanh[/?]) 



(60) 



327r/3 2 s 2 m 2 

where is defined above. Note that if more than one KK neutrino species is present, 
two related processes will also take two neutrinos of different species into their two zero 
modes, or one neutrino and one anti-neutrino of different flavors into their zero modes. 
In both cases, we have a single t-channel Feynman diagram and the cross sections are, 

gi(4s - 3m 2 ) 



V\V2) 



647r/3sm 2 ' 
gj ((3(4s + 9m 2 ) + 8m 2 L) 
64irf3 2 sm 2 



(61) 
(62) 
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Figure 11: Feynman diagrams for annihilation into two zero mode leptons vv. 

Initial state Final state Feynman diagrams 

(zero modes of SM fields) 



J 1 ) F M 

e R e R 



q q s(B<®) 
v V s(B^) 

e~ e+ s(-B<°>), 
(f) 0* s(B^) 
Z Z f( e g>), u{ef) 

7 7 t(eg>), u(eg>) 

Z 7 t(eg>), «(eg } ) 

eg } eg } e" e" t(# (1) ), w(S (1 )) 

eg } BO) e~ 7 S (e~), t(eg } ) 

e" Z s(e-), t(eg>) 



Table 3: Same as Table p] but for coannihilation of e 



(i) 

R ■ 



with = g z . Finally, we can have cross-flavor transition between into charged 

lepton zero modes, e^e^ through a t-channel W± exchange. The corresponding cross 
section is given by the result for a(u^ i4 ~~ * ^1^2) with g L = e/y/2sw 

C Coannihilation Cross Sections 

The eg can annihilate with eg into quark (and other family lepton) zero modes through 
an s-channel B^ ' (Figure^). The cross section is given by, 

m m N c gfY e 2 ( Yf + Y 2 )(s + 2 m 2 ) 
o&M^ff) = R[ 2A jj (63) 

Note that since eg^ is a weak singlet, this formula also applies for annihilation into zero 
modes of the neutral lepton partner. 
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Annihilation into e + e zero modes proceeds through an s-channel B^ or by exchang- 
ing a t-channel (Figure [7]). The cross section is, 

(1) _ (1) + gft%[5fo + 2(2 3 + 3m 2 )L] gft? a [g (4. + 9m 2 ) + 8m 2 L] 

^ 6 6 j " 32vr / 3 2 S 2 + 64vrm 2 /? 2 S 

+ (^ + 2m 2 ) 
+ 24^ (64) 

where, L is defined in Section [B|. Note that we have chosen to include the decay into 
left-handed electrons in this result, though we could have equally well considered it part 
of Eq. (p3"|). Evident from comparison of the two equations, this is simply a matter of 
book-keeping. 

Annihilation into zero modes of Higgs bosons occurs through an s-channel B^ (Fig- 
ure |8|) with cross section, 

where the factor of 2 to sum over the two entries of the doublet is included. 

Annihilation into ZZ, Zj, and 77 are mediated by t- and w-channel (Figure |9[). 
This process can be more conveniently described as the single channel e^e^ — > B^B^°\ 
equivalent to the sum of these three 7 and Z processes. The cross section is, 



ale 



(i)_(i) R ( 0)R (oh 9$K ( 2 [ g2 + 4 ^ 2g ~ 8 ^ 4 ] ArcTanh[/3] - (3s\s + 4m 2 ]) 
R e R -> B B ) - 



Annihilation into W + W zero modes is zero in the limit in which one ignores EWSB 
effects, because is a singlet. 

Finally we have the process e^e^ — > e~e~ ( eg' ) eg' ) — > e + e + ) via exchange of a t- or 
it-channel B^\ This cross section is 

{1) (1) _ ^n 4 fl (^(23 - m 2 ) + 2m 2 (4s - 5m 2 )ArcTanh[/3]) 
a(e « e ^ e e ) = ' (67) 

for lepton KK modes of the same flavor, and, 

f W W - ^4(4* -3m 2 ) 

647rpsm z 

/ (D-d) - -k ^4 + 9m 2 ) + 8m 2 L) 

o-(ek Vk e M ) = ~ ^7-^ 2 69 

647rp^sm^ 

for two modes of different lepton flavor. 

Coannihilation of a B^ with e R (or e^) proceeds into either e~ Z or e~ 7. In the 
limit in which the Z mass is disregarded, we can equally well describe this as a single 
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Figure 12: Feynman diagrams for B^f^ annihilation into a zero mode / and vector 
boson. 



channel into eB^. The cross section is given by, 



a(B^ 



B<®e-) = 



Y e 4 R gf[(3s(2m 2 - s) + m 2 (6m 2 - 2s) L] 
48irf3 2 s 2 m 2 



9Q7r/3sm 2 

Y e R 9t \Ps(s - Am 2 ) - 2m 2 (s + 6m 2 ) L] 
967i/3 2 s 2 m 2 



(70) 



The KK modes of the left-handed electron are somewhat more complicated because 
they involve the SU(2) bosons as well as the U(l) boson. We find it convenient to consider 
neutral zero mode gauge bosons in the W^°\ B^ basis. Annihilation into fermions is given 
by 



N c g i (s + 2m 2 ) 
24tt(3s 2 



where the coupling, 



g 2 = (gfY f Y eL 



9l 



(71) 



(72) 



is in terms of the hypercharges (Y) and third component of weak iso-spin (T 3 ) for the 
fermion and the left-handed electron. For annihilation into zero modes of electrons of the 
same family we also have t-channel exchange of -B*- 1 -* and W^°\ and for annihilation into 
zero modes of the neutral lepton partner we have t-channel W±^ exchange. Both results 
may be expressed, 



e+e _ ^ g 2 L g 2 [5(3s + 2(2s + 3m 2 )L] + g\ [(3 (4s + 9m 2 ) + 8m 2 L] 



32tc(3 2 s 2 
g A (s + 2m 2 ) 
24irf3s 2 



64 7r m 2 j3 2 s 



(73) 



where g 2 is defined above, for the specific case of / an electron or neutrino, and g^ 
e/2s\yCw f° r electrons and = e/^/2s w for the neutrino final state. 
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Initial state Final state Feynman diagrams 

(zero modes of SM fields) 



(1) 




q q 
v V 

e~ e + 
<p <p 

z z 


S (Z, 7 ),t(£« W 3 (1) ] 

S 17> Z J 














7 7 






lZ 


t(ei 1} ), u(e£>) 






W + W~ 


s( 7 ,Z),«(i/W) 


JXi 

e L 




e 6 


+CR(1) M/ (1 h 


e L 




6 /i 


+ ( 'r(i) roWl 


$ 


e (1) 


/x _ e + 


l \ n 1 vv ± ) 


J 1 ) 
e L 


77(1) 

IS 


y y 

z w- 

7 W~ 


s(W ) 
s(W-),t(B^\ W ^) 

s(WJ) 
s(W_), t(e£>), u{vX) 
s(W-), t(e ( l ] ), u{vM 








t(BW,wP), u(W m 




v w 




u(W { ± ] ) 






/i _ 77 


t(BW, W P) 



Table 4: Same as Table p] but for coannihilation of . 

Under our approximation in which EWSB effects are neglected, annihilation of 
into W bosons is equal to the process v^V^ — > W + W~ given in Equation (|57"D , by SU(2) 
invariance. Similarly, the sum of the annihilation processes into ZZ, 7 Z, and 77 are also 
equal to the process pX)jyW — > ZZ given in Equation (|56|) and the sum of annihilation 
into both components of the Higgs doublet is given by the process i/^'pi 1 ) — > (jxp* in 
Equation (54]). Furthermore, like-sign annihilation e£ e~e~ is equal to 



vv, and thus is given by Equation (BO). 

e L annihilate into zero modes of quarks or other family leptons through an s- 
channel with cross section given by Equation ([H]), replacing with the appropriate 
charged current coupling: g 2 — > g\j1. Annihilation into leptons of the same family also 
includes t-channel exchange of Z^\ The cross section is given by Equation (|73| ) with 
g 2 — > g\j2,g\ — > (e/2swCw) 2 [~ 1/2 + Syp]- Annihilation into Higgs occurs through an 
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s-channel and may be obtained from Equation (|54j) with the replacement g\g\ — > 
e 4 /4sy V . Annihilation into gauge bosons is most simply described in terms of annihilation 
into W-B<®, and W~W^ . The first process is mediated by i-channel e\} exchange and 
it-channel exchange, and is given by Equation ( |56| ) with g\ — > — e 2 /2\/2swCw- The 
second is simply obtained from — > W + W in Equation (|57f) by SU(2) invariance. 
Finally, e^u^ 1 — > e^z/ is mediated by t-channel Z^ 1 ) and w-channel W^/ exchange, and 
is given by, 

(i) (1) _ g 2 gl [4m 2 (4s - 5m 2 )ArcTanh[/3] + m 2 (3s] 

a(e L v - e i/) _ 

+ fo[(g t 4 + ^)(4 a -3m 2 )] 
647r/3 2 s 2 m 2 

with s 2 = e 2 /(2s 2 / c 2 y )(-l/2 + s 2 w ) and <? 2 = e 2 /2s 2 / 

If there are multiple families of KK left-handed electrons or neutrinos, there will also 
be flavor-changing annihilation of p^v^ — > jt/^i/ 1 ) — > z-^e - , and /i^ X^ 1 ) — > /i~z7. 
The cross sections for the first two processes are both given by Equation (|6T|), with the 
replacements g\ — > e 2 /(2s 2 v c 2 v )(—l/2 + s 2 ^), and g\ — > e 2 /2s 2 v . The cross section for 
/i^ X^ 1 ) — ► n~V is given by Equation (|62[) with the replacement g 2 — * e 2 /(2s 2 v c 2 v )(— 1/2 + 
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